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ABSTRACT 

We present X-ray, optical, near-infrared (IR), and radio observations of GRBs 110709B and 111215A, as 
well as optical and near-IR observations of their host galaxies. The combination of X-ray detections and deep 
optical/near-IR limits establish both bursts as "dark". Sub-arcsecond positions enabled by radio detections 
lead to robust host galaxy associations, with optical detections that indicate z < 4 (110709B) and z « 1.8-2.7 
(11 1215A). We therefore conclude that both bursts are dark due to substantial rest-frame extinction. Using the 
radio and X-ray data for each burst we find that ORB 1 10709B requires A^"""^ > 5.3 mag and GRB 1 1 1215A 
requires Ajj?^' ^8.5 mag {z = 2). These are among the largest extinction values inferred for dark bursts to 
date. The two bursts also exhibit large neutral hydrogen column densities of A^n.int ^ 10^^ cm~^ (z = 2) as 
inferred from their X-ray spectra, in agreement with the trend for dark GRBs. Moreover, the inferred values 
are in agreement with the Galactic Ay -Nh relation, unlike the bulk of the GRB population. Finally, we find 
that for both bursts the afterglow emission is best explained by a collimated outflow with a total beaming- 
corrected energy of E-^ + Ek ~ {1 - 9) 10^' erg (z = 2) expanding into a wind medium with a high density, 
M ss (6-20) X 10"^ M© yr"' (n « 100-350 cm""* at « 10'^ cm). While the energy release is typical of long 
GRBs, the inferred density may be indicative of larger mass loss rates for GRB progenitors in dusty (and hence 
metal rich) environments. This study establishes the critical role of radio observations in demonstrating the 
origin and properties of dark GRBs. Observations with the JVLA and ALMA will provide a sample with sub- 
arcsecond positions and robust host associations that will help to shed Ught on obscured star formation and the 
role of metallicity in GRB progenitors. 

Subject headings: gamma rays: bursts — dust — extinction 



1. INTRODUCTION 

Long-duration gamma-ray bursts (GRBs) have been linked 
to the deaths of massive stars, and hence to star forma- 
tion activity, through their associatio n with star forming 
galaxies (e.g., fPiorgovski et al. 1998; ' Fruchter et al.l l2006t 
IWainwright et al. 2007) and with Type Ic supernova explo- 
sions (e.g., Woosle y & Bloom 2006). Across a wide range 
of cosmic history a substantial fraction of the star forma- 
tion activity (~ 70% at the peak of the star formation history, 
z ~ 2-4) is obscured by dust, with about 15% of the total star 
formation rate density occurring in ultra-luminous infrared 
galaxies (e.g Bouwens et al. 2009; Reddy & Steidel 20091 
iMurphv et al.ll201 iT) . As a result, we expect some GRBs to 
occur in dusty environments that will diminish or completely 
extinguish their optical (and perhaps even near-IR) afterglow 
emission. Such events can be used as signposts for the loca- 
tions and relative fraction of obscured star formation across a 
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wide redshift range. In addition, they can provide insight into 
the role of metallicity in GRB progenitors since dusty envi- 
ronments generally require substantial metallicity. 
These so-call ed optically - dark GRBs are indeed known 



to exist (e . g.. lGrootetal.1 1 19981; iDiorgovski et al.1 120011; 



(e . g.. luroc 
^ I200U iPi] 



> et al.1 



20021) . but the lack of an opti- 
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cal detection does not necessarily point to dust obscuration. 
Most prosaically, the lack of detected optical emission may 
be due to inefficient follow-up observat ions, or to intrinsi- 
cally dim events (e.g.. iBerger et al.ll2002h . Another potential 
origin of dark bursts is a high redshift, with the optical emis- 
sion suppress ed by Lya absorption at Apbs 5, 1216A x(l + z) 
jHaislip et a l. 2006; Salvaterra et al. "2009; " ^vir et al.ll2009l: 
ICucchiaraet al. 201 li) . Such events are clearly of great inter- 
est since a dropout above ^ 1 jim points to z > 7, while events 
which are also dark in the near-IR can potentially arise at 
z > 18 (/T-band dropout). Naturally, high redshift bursts will 
also lack host galaxy detections in the optical band. On the 
other hand, their afterglow emission redward of the Lya break 
will follow the expected synchr otron spectrum (/v k 
with /3 w 0.5 - 1 (ISari et al.lll998l) . 

To determine whether a burst is genuinely dark, due to 
extinction or to a high redshift rather than to an inefficient 
search, it is common to compare the observed limits with the 
expected optical/near-IR brightness based on the X-ray after- 
glow. This approach relies on the simple power law shape of 
the afterglow synchrotron emission. One definition of dark 
bursts uses an optical to X-ray spectral index of /3ox J; 0.5 
(/3 is the slope of the synchrotron spectrum in log space; see 
iJakobsson et al.l [20041 for details), since this is the shallow- 
est expected slope in the standard afterglow model. A vari- 
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ant of this condition uses knowledge of the X-ray spectral 
index (0x), and defines bursts as dark if /?ox -fix ^-0-5 
jvan der Horst et al] 12009 '). since this is the shallowest ex- 
pected relative slope. These definitions can reveal evidence 
for dust extinction even if an optical afterglow is detected. 

This approach has been used by several gr oups to identify 
and stu dy dark bursts. In the pre-Swift sample lJakobsson et alj 
(120041) fou n d that about 10% of the bursts are dark, while 
iKann et all (120061) found a mean extinction o f A\°'^ » 0.2 
mag, indicating that most events are not dark. ISchadv et al.l 
12007) studied several bursts with Swift X-ray and UV/optical 
detections and found a mean extinction level of Ay°^' w 0.3 
mag, and that bursts which are dark blueward of V-band are 
likely to ha ve rest-frame extinc tion about an order of magni- 
tude larger. iMelandri et al.l (120081) used rapid optical observa- 
tions of 63 Swift bursts and found that about 50% showed ev- 
ide nce of mild extincti on. A similar conclusion was reached 
by ICenko et al.l ( l2009 l) based on rapid optical observations 
of 29 Swift bursts, with an 80% detection fraction, but with 
roughly half exhibiting suppression with respect to the X- 
ray emission. A f ollow-up study of the latter sample by 
iPerlev et"an (l2009ah aimed at identifying host galaxies in the 
optical (thereby ruling out a high redshift origin) and found 
that < 7% of Swift bursts are located at z > 7. The major- 
ity of the dark bursts in their sample instead require Aj y^' > 1 
mag, with a few cases reaching ~ 2-6 mag. Melandr i et al.l 
(j2012) studied a complete sample of 58 bright Swift bursts, 
of which 52 have known redshifts and found that the fraction 
of dark bursts is about 30%, mainly due to extinction (nearly 
all the dark bursts in their sample have z < 4). Thus, the af- 
terglow emission of ^ 1/3 - 1/2 of all GRBs are affected by 
dust, although in mos t cases the required extinction is mod- 
est, A^r' ~ 0.3 mag (Kam et all 120061 ISchadv et al.l 120071: 
iPerley et"ar. 2009b; Greiner et al. 2011). There are only a few 
known cases with large extinction of A^°^' ^ few mag. 

Concurrent studies of the neutral hydrogen column den- 
sity distribution, inferred from the X-ray afterglow spectra, 
suggest that dark bursts exhibit systematically larger values 
of A^H.int }t 10^^ cm~^ compared to bursts with little or no 
extinction, which have a med ian of A^H.int ^ 4 x 10^' cm~^ 
dCampana et al. 2012; Margut tietan 120121) . On the other 
hand, the measured extinction for GRBs is generally lower 
than expected based on the values of A^H,int and the Galactic 
Ay—NH relation, previously attribut ed to dust destruction by 
the bright X-ray/UV emission (e.g.. Ga lama & Wiiersll20()l1; 
Hchady et al. 2007; Perley et al. 26o9b). 

Finally, the host galaxies of at least some dark bursts appear 
to be redder, more luminous, more massive, an d more metal 
rich than the hosts of optically-bright GRBs dBerger et al. 
2007"; 'Perlev et al."2009b( iLevesque et ani2010l;1Perlev et al. 
201 la; Kriihler et al. 201 f), potentially indicating that the ex- 
tinction is interstellar in origin rather than directly associ- 
ated with the burst environment. One specific example is 
GRB 0301 15, not only one of the first examples of a heavily 
extinguished GRB afterglow, but residing in a host classified 
as an extremely red object (ERO; Levan et al. 2006). Some 
host galaxies of optically-bright bursts have been detected in 
the radio and (sub)millimeter ranges (e. g.. .Tanvir et al.ll2004i 
indicat ing obscured star f ormation rates of ~ 10^-10^ Mq 
yr"' ( Berg er et al.l l2003ah . This suggests that the dust dis- 
tribution within the host galaxies is patchy and that not all 
GRBs in dusty host galaxies are necessarily dust-obscured 
dBerger etalj|2003at iPerlev et al.ir2009b,) . Further work by 



ISvensson et al.l d2012l) indicates that dark GRB host galax- 
ies may be systematically redder and more massive than op- 
tically bright GRB host galaxies and have lower metallicity 
than sub-millimeter galaxies, leading to prior biases in GRB 
host galaxy population studies. 

Here we present multi-wavelength observations that reveal 
two of the darkest known bursts to date, GRBs 110709B 
and 111215A. These events have such large rest-frame ex- 
tinction that they lack any afterglow detection in the optical 
and near-IR bands (to Aobs ~ 2.2 fim) despite rapid follow- 
up. However, both events are detected in the radio, pro- 
viding sub-arcsecond positions and secure associations with 
host galaxies detected in the optical. This allows us to rule 
out a high redshift origin, while the combination of radio 
and X-ray data place robust lower limits on the rest-frame 
extinction, with Ay°'*' > 5.3 and > 8.5 mag (for z = 2), re- 
spectively. The radio and X-ray data also allow us to deter- 
mine the explosion properties and the circumburst densities. 
This study demonstrates that radio detections with the JVLA, 
and soon the Atacama Large Millimeter/submillimeter Array 
(ALMA), provide a promising path to accurate localization of 
dark GRBs, robust estimates of the extinction in the absence 
of any optical/near-IR detections, and a comparative study of 
their explosion properties and parsec-scale environments. 

The plan of the paper is as follows. We describe the multi- 
wavelength observations of GRBs 1 10709B and 1 1 12 15 A and 
their host galaxies in ^ These observations establish that 
both events are dark due to extinction. We model the ra- 
dio and X-ray data to extract the minimum required extinc- 
tion, as well as the explosion and circumburst properties in 
^ In ^ we place the large inferred extinction values and 
the inferred burst properties within the framework of existing 
samples of optically-bright and dark bursts. We summarize 
the results and note key future directions in ^ Through- 
out the paper we report magnitudes in the AB system (un- 
less otherwise noted) and use Galactic extinction values of 
E(B-y) « 0.044 mag for G RB 1I0709B and E(B-V) » 
0.057 mag for GRB 111215A dSchlaflv & FinkbeineiillOlTI) . 
We use the standard cosmological parameters: Ho = 7l km s~^ 
Mpc-i, flA = 0.73, and ilu = 0.27. 

2. OBSERVATIONS 
2. 1 . Discovery and Burst Properties 
2.1.1. GRB110709B 

G I^ I10709B was discov ered on 2011 July 9 at 21:32:39 
UT (Cumm ingsetalJl20Ilal) by the Swift Burst A lert Tele- 
scope (BAT; 15-150 keV: Barthelmv et al J l2005h, and b y 
Konus-WIND (20 keV-5 MeV; Golenets kfietal] l20Tlh 
Swift X-ray Telescope (XRT; 0.3-10 keV: 'Burrow s et alJ 
2005) observations commenced at 5t JCumming s et alJ 

201 lai) . locaHzing the X-ray afterglow to RA=10''58'"37.08^ 
Dec=-23°27'17.6" (J2000), with an uncertainty of 1 .4" (90% 
containment, UVOT-enhanced; Beardmore et al. 2011a). 
Swift UV/Optical Telescope (UVOT: .Roming et al.i.200 5) ob- 
servations began at 5t 90 s, but no afterglow candidate 
was identifi ed to a 3g h mit of « 2 1 . 1 mag in the white filter 
( |Cu mming s et alJl201 lah . A large flare was detected with the 
BAT, XRT, and Konus-WIND at (5f fti 1 1 min with an in tensity 
c omparable to the main event ( iBarthelmv et alj|201 lah . 

IZhang et"an (^012) presents a detailed analysis of the 
prompt emission. The initial trigger had a duration and flu- 
ence in the BAT 15-150 keV band of Tgo « 56 s and = 
(8.95!|]:^^) X lO"*" erg cm'^, while the flare had Tgo « 259 s 
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and = (1.34!° °^) x 10"^ erg cm-^ (I Zhang et alJlMTl) . No 
emission was detected in the BAT band at 6t ~ 180-485 s. 
Combining both events, the resuhing duration is Tgo « 846 ±6 
s (Cummi ngs et al.ll201 Ibl) . The initial trigger had a fluence 
in the Ko nus-WIND20 -5000 keV band of (2.6 ±0.2) x 10"^ 
erg cm"^. IZhang et"an (12012) calculate the total fluence in the 
SwiftfXRT band to be 4.07 ± 0.56 x 10"'' erg cm'^. 

2.1.2. GRB111215A 

GRB111215A was discove red on 2011 Dece mber 15 at 
14:04:08 UT by the 5w/^/BAT (lOates et al.l2011h . SwiftlXKV 
observations began at 5t w 409 s and localized the X-ray 
afterglow to RA=23''18™13.29\ Dec=+32°29'38.4", with 
an uncertainty of 1 A" ( 90% containment, UVOT-enhanced; 
iBeardmore et al.ir201 Ibl) . No afterglow candidate was iden- 
tified by th e UVOT to a 3cr limit of > 21.7 mag in the 
white filter (IOatesll2()Tl1) . The duration and fluence are Tgo « 
796 s and = (4.5 ± 0.5) x 10"^ erg cm'^ (15-150 keV; 
lBartheknvet'ani2011bh . 

2.2. X-ray Observations 
2.2.1. GRB110709B 

We analyzed the XRT data using the HEASOFT package 
(v6.11) and corresponding calibration files. We utilized stan- 
dard filtering and screening criteria, and gen erated a count- 
rate hg ht curve following the prescriptions bv iMargutti et alj 
^W). The data were re-binned with the requirement of a 
minimum signal-to-noise ratio of 4 in each temporal bin. 

The XRT data exhibit significant spectral evolution during 
the early bright phase (5t < 1 ks), with spectral hardening 
when the burst re-brightened at ^ 1 1 min. Due to this varia- 
tion, we performed a time-resolved spectral analysis, accumu- 
lating signal over time intervals defined to contain a minimum 
of about 1000 photons. Spectral fitting was done with Xspec 
( V 1 2 . 6 : . Arnaud 1996) . assuming a photoelectrically absorbed 
power law model and a Galactic neutral hydrogen co lumn 
density of A^h.mw = 5.6 x lO^" cm'^ jKalberla et alj2005h . We 
extracted a spectrum in a time interval when no spectral evo- 
lution was apparent (15-150 ks) and used this to estimate the 
contribution of additional absorption. We find that the data 
are best modeled by an absorbed power-law'' model with a 
spectral photon index of F = 2.3 ±0.1 and excess absorption 
of A^H.int = (1.4 ±0.2) X 10^' cm-2 at z = (90% c.l., C-stat 
= 559 for 604 degrees of freedom). The excess absorption 
was used as a fixed parameter in the time-resolved spectral 
analysis, allowing us to derive a count-to-flux conversion fac- 
tor for each spectrum. The resulting unabsorbed 0.3 - 10 keV 
flux light curve properly accounts for spectral evolution of the 
source. We note that for the purpose of display in Figure [T] 
we binned the XRT data by orbit at (5f < 2 d and by multiple 
orbits thereafter 

We also observed ORB 110709B for 15 ks with the Ad- 
vanced CCD Imaging Spectrometer (ACIS-S) on-board the 
Chandra X-ray Observatory on 2011 July 23.60 UT {5t « 
13.7 d). The afterglow is detected with a count rate of 
(4.1 ±0.5) X 10"^ s"' (0.5-8 keV), coiTesponding to an un- 
absorbed 0.3-10 keV flux of (6.0 ± 0.7) x lO"'"* erg s"' cm'^ 
using the XRT spectral parameters. Relative astrometry with 
the Chandra position will be discussed in §2.6. A second, 9 
ks Chandra/ ACIS-S observation was performed on 201 1 Oc- 
tober 31.83 UT (^f « 113.4 d). The X-ray afterglow is not 

' ' tbabs x ztbahs X pow 



detected to a 3(7 limit of < 9 x 10"'^ erg s"' cm"^, well above 
the extrapolation of the early light curve to this epoch (Fig- 
ure [T]!. 

At (5f > 12 d the XRT light curve flattens to a level of about 
lO"'"* erg s"' cm"^. The Chandra observations reveal that this 
is due to a contaminating source located about 3.8" from the 
GRB position, well within the point spread function (PSF) of 
XRT. We subtract the flux of this source, (6.0 ±0.7) x lO"'"* 
erg s"' cm-2, from the XRT light curve of GRB 1 10709B in 
Figure [T]'^ 

2.2.2. GRB 1 1121 5A 

We analyzed the XRT data for GRB 1 1 12 15 A in the same 
manner described above. At St < 2 ks the light curve ex- 
hibits flares with a hard-to-soft spectral evolution. To ob- 
tain a reliable estimate of the intrinsic neutral hydrogen ab- 
sorption in addition to the Galactic value (A^h.mw = 5.5 x 10^° 
cm -2; iKalberla et al]l2005l) . we extracted a spectrum collect- 
ing the photon count data in the time interval 5 - 2000 ks, 
when no spectral evolution is apparent. The data are best 
modeled by an absorbed power-law with T = 2.2 ±0.1 and 
A^H.int = (3.1 ±0.4) X lO^i cm-2 at z = (90% c.l., C-stat = 
499. 13 for 512 degrees of freedom). The resulting unabsorbed 
0.3 - 10 keV flux light curve is shown in Figure |2l 

2.3. Optical/Near-IR Afterglow Limits 
2.3.1. GRB110709B 

We observed GRB 1 1 0709B with t he Ge mini Multi-Object 
Spectrograph (GMOS; iHook et alj |2004 on the Gemini- 
South 8-m telescope in r-band on 2011 July 10.03 and 13.96 
UT (St w 3.2 hr and w 4.1 d). We analyzed the data with 
the gemini package in IRAF, but did not detect any sources 
within the XRT error ckcle to 3a limits of > 23.8 mag 
(5t « 3.2 hr) and > 25.1 mag (St « 4.1 d). A comparison of 
the early r-band limit (/v.opt ^1-1 mJy) to the X-ray flux den- 
sity at the same time (iv x ~ 4.0 /iJy) indicates an optical to 
X-ray spectral index /3ox ^ -0.2. This is substantially flatter 
than the minimum expected spectral index of /3ox = 0.5, in- 
dicating that GRB 1 10709B is a dark burst; for /3ox = 0.5 the 
expected r-band brightness is « 19 mag. Similarly, relative 
to the X-ray spectral index we find /3ox- (3x ^ -1-5, clearly 
satisfying the dark burst condition (iox- ^ -0.5. 

We also observed the burst with the Gamma-Ray Burst 
Optical/Near-lnfrared Detector (GROND: lGreiner et alJ2"007l 
2008) mounted on the Max Planck Gesellschaft / European 
Southern Observatory 2.2-m telescope at La Silla Observa- 
tory starting on 201 1 July 9.94 UT {5t « 1 hr) simultaneously 
in g'r'i'z'JHK^, with an average seeing of 1.5" and an av- 
erage airmass of 1.3. We analyzed the data using standard 
pyraf /IRA F tasks, following the procedures described in 
iKriihler et alJ ([2008). We calibrated the optical channels with 
zero-points computed during photometric conditions, and the 
near-lR channels using stars in the 2MASS catalog. Calibra- 
tion uncertainties vary in the range 0.02-0.12 mag for JHKs, 
and we add these systematic errors with the statistical errors 
in quadrature. We do not detect any source in the individual or 
stacked images within the XRT error circle. The upper limits 
are listed in Table [T] derived by forcing the photometry at the 
position of the radio counterpart 02.4b . A comparison of the 

We note that the flux from the contaminating source is constant between 
the two epochs of Chandra, and rules out any origin as a lensed counterpart 
of the GRB, which may have caused the two 7-ray triggers. 
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Ks band limit at (5f w 2.7 hr (/v.nir ^ 69 /iJy) to the X-ray flux 
density at the same time indicates a near-IR to X-ray slope of 
< 0.35, indicating that GRB 110709B is a dai'k burst in the 
near-IR as well. 

The absence of optical/near-IR emission can be due to dust 
extinction, or alternatively to a high redshift, z > 18 based on 
the /Tj-band non-detection. We rule out a high redshift origin 
ing22] 

2.3.2. GRB1U215A 

Several optical afterglow searches on timescales of 8t ~1 
min to 6.8 hr led to non-detections of an afte rglow with 
lim its of > 17 5-22.8 mag in va-ious filters ('Xm '2011'; 
XuHIOll lUsuiljgOn I^nd"evll201 U iGorbovskov e t al. 2011; 



Aceituno et al.ll201ll : lRumyant sev 201 1). The deepest opt ical 
limit is iriR > 22.8 mag at a mid-time of St « 37 min dXul 
1201 lb . corresponding to a flux density of /V opt 5, 2.3 nJy. 
Near-IR observations al so led to non-detections (iD^Avanzol 
1201 It iTanvir et"ani201 Ih . with the deepest liniit bein g m/j- > 
19.3 mag (Vega) at St ^ 0.20 d jTanvir et al.ll2011h . coitc- 
sponding to /v,nir ^12 fiJy. 

A comparison of the deepest optical limit to the X-ray flux 
density at the same time (/v.x ~ 54 /iJy) indicates /3ox ^ -0.5. 
Similarly, a comparison of the deepest near-IR limit to the 
X-ray flux density at the same time (F^ x ~ 4.2 /iJy) indi- 
cates /3nirx ^ 0.15. In comparison to the X-ray spectral in- 
dex we find Pox- ^ and /3nirx-/3x ^ ~1-05. Thus, 
GRB 111215A is a dark burst in the optical and near-IR. As 
in the case of GRB 110709B, the absence of optical/near-IR 
emission can be due to dust extinction, or alternatively to a 
high redshift, z > 18 based on the Aj-band non-detection. 

2.4. JVLA Centimeter Observations 
2.4.1. GRB110709B 

We observed GRB 1 10709B with the NRAO Karl G. Jansky 
Array (JVLA) beginning on 201 1 July 1 1.97 UT {St « 2.1 d) 
at a mean frequency of 5.8 GHz and detected a single, unre- 
solved radio source within the XRT error circle. This detec- 
tion provided the first accurate position for the burst. Follow- 
up observations demonstrated that the source initially bright- 
ened and subsequently faded away, establishing it as the radio 
afterglow of GRB 1 10709B. 

All observations utilized the WIDAR correlator 
(iPerley et al.l 1201 Ibl) with - 2 GHz bandwidth. We cal- 
ibrated and analyzed the data using standard procedures jn 
the A stronomical Image Processing System (AIPS; iGreisenI 
l2003h . We excised edge channels and channels affected by 
radio frequency interference, reducing the effective band- 
width by ^ 25% at 5.8 GHz. Due to the low declination of the 
source we also excised data when > 2 m of a given antenna 
was shadowed by another antenna. For the observation at 
21.8 GHz we performed reference pointing at 8.4 GHz and 
applied the pointing solutions, per standard high frequency 
observing procedures. We observed 3C286 for band-pass and 
flux calibration and interleaved observations of J 1048- 1909 
for gain calibration every 3 m at 21.8 GHz and Jl 112-2158 
every 4 m at 5.8 GHz. The resulting flux densities at 5.8 
and 21.8 GHz are listed in Table |2] The uncertainties are la 
statistical errors, and we note an additional uncertainty in the 
absolute flux scaling of ^ 5%. The 5.8 GHz light curve is 
shown in Figure [T] 

Finally, to determine the position of the radio afterglow 
we fit the source in each image with a Gaussian profile 



(AIPS task JMFIT), and calculated the mean radio position 
weighted by the resulting statistical uncertainties. We find 
a weighted mean position of RA=10*'58™37.113'* (±0.001), 
Dec=-23°27'16.76" (±0.02). 

2.4.2. GRB 11 121 5A 

We observed GRB 1 1 1215A with the JVLA beginning on 
2011 December 17.00 UT (St w 1.4 d), and subsequently de- 
tected the radio afterglow starting at (5f w 3.35 d. Observations 
between 1.4 and 140 days were obtained at mean frequencies 
of 5.8, 8.4 and 21.8 GHz. We used 3C48 for band-pass and 
flux calibration and interleaved observations of J23 11+3425 
every ^4 m for gain calibration at 21.8 GHz and J2340+2641 
every ^5.5 m for gain calibration at 5.8 and 8.4 GHz. The 
data were analyzed in the same manner described above, and 
the resulting flux densities are listed in Table [3] The light 
curves are shown in Figure s|2] and [3] 

We obtain a weighted mean position from all 21.8 
GHz detections of RA=23*'18™13.314-^ (±0.004), 
Dec=+32°29'39.07" (±0.06). The JVLA data (along with 
CARMA observations: ^2.5b provide the only sub-arcsecond 
position for GRB 1 1 12 15 A. 

2.5. Millimeter Observations 

2.5.1. GRB110709B 

GRB110709B was observed at a frequency of 345 
GHz with the Large APEX Bolometer Camera (LABOCA; 
Siring o et al.l 1200 9) on the Atacama Pathfinder Experiment 
(APEX) telescope beginning on 201 1 July 1 1 .90 UT (St « 2.0 
d). No submillimete r counterpart was detected to a 3a upper 
limit of < 6.9 mJy (ide Ugarte Postigo et all201 Ih . 

2.5.2. GRB 11 121 5A 

We o bserved GRB 1 1 12 15 A with the Submiflimeter Array 
(SMA; iHo et"an 12004 ) beginning on 2011 December 18.15 
UT (St « 2.7 d) at 230 GHz with 8 GHz bandwidth. We 
performed gain calibration using J231 1+344 and J2236+284, 
band-pass calibration using 3C84, and flux calibration using 
real-time measurements of the system temperatures, with ob- 
servations of Uranus and Callisto to set the overall flux scale 
(accurate to about 10%). We calibrated the data using the 
MIR software package, and analyzed and imaged the obser- 
vation with MIRIAD. We do not detect the GRB afterglow to 
a 3a limit of < 2.6 mJy. 

We also observed GRB 11 12 15 A with the Combined Ar- 
ray for Research in Millimeter Astronomy (CARMA) begin- 
ning on 2011 December 16.97 UT (St « 1.38 d) at 93 GHz, 
and continued observations for « 16 d; see Table [3] We 
used Neptune as the primary flux calibrator, and J 1824+568 
and J1638+573 as band-pass and gain calibrators, respec- 
tively. The overall uncertainty in the absolute flux calibra- 
tion is ^ 15%. We calibrated and visuahzed the data with 
the MIRIAD software package (ISaultet al.lll995Ll201 Ih . and 
measured the flux densities with AIPS for consistency with 
the JVLA observations (TableO. 

We fit the position of the detected mm afterglow in each 
epoch with a Gaussian fixed to the beam size, and com- 
pute the mean position, weighted by the uncertainty in 
each epoch, leading to RA=23^18'"13.32P (±0.007) and 
Dec=+32°29'39.04" (±0.07). This position is consistent wifli 
the Swift/X-iay position. 

2.6. Host Galaxy Observations and Redshift Constraints 
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2.6.1. GRB110709B 

We obtained Hubble Space Telescope (HST) observations 
of GRB110709B with the Wide Field Camera 3 (WFC3; 
iKimble et al.ll2008l) on 2011 November 8.94 UT (6t « 122 
d) using the ultra-violet imaging spectrograph (UVIS) with 
the F606W filter, and on 201 1 November 12.94 UT (St « 126 
d) using the infrared (IR) channel with the F160W filter (PI: 
Levan). One orbit of observations was obtained in each fil- 
ter, with exposure times of 2610 s (F606W) and 2480 s 
F160W). The data were processed using mult idrizzle 
Fruchter & Hookll2002l: iKoekemoer et alj|2006h with output 
pixel scales set at 0.02" and 0.07", respectively. 

To locate the absolute radio and X-ray afterglow positions 
on the HST images we tie the F606W and F160W images 
to the 2MASS reference frame using a wider field Gemini- 
South GMOS r-band image as an intermediary. Using 20 
common sources with 2MASS, the absolute astrometry of 
the Gemini image has an rms scatter of (Tgmos-2mass =0.15" 
in each coordinate. The HST images are tied to the GMOS 
image using 30 common sources with a resulting rms scat- 
ter of (Jhst-gmos = 0.03" in each coordinate. We further re- 
fine the Chandra afterglow position relative to the HST as- 
tro me trie frame using a single common source. This leads to 
a shift in the X-ray position of SRA = -0.10" (-0.07") and 
(5Dec = +0.20" (+0.18") relative to the F606W (F160W) ref- 
erence frame. There are no common sources between the HST 
and JVLA images. 

The resulting uncertainty in the radio afterglow position in 
the //5r reference frame is dominated by the absolute astrom- 
etry, and corresponds to a radius of 0.22" (la; the centroid 
uncertainty of the radio position is only 0.02" in each coor- 
dinate). The uncertainty in the X-ray afterglow position is 
reduced by the relative tie of the Chandra and HST reference 
frames, leading to a radius of 0.10" (la; the centroid uncer- 
tainty of the Chandra position is about 0.08" in each coordi- 
nate). Within the uncertainty regions we find a single galaxy 
in both //5r images (Figure |4]i, which we identify as the host 
of GRB 110709B. Photometry of this source using the tabu- 
lated WFC3 zeropoints gives mpeoew = 26.78 ±0.17 AB mag 
and mpieow = 25 . 1 3 ± 0. 10 mag, corrected for Galactic extinc- 
tion. 

To determine the probability of chance coincidence, fcc(< 
SR) = 1 - exp[-7r((5/?)^cr(< m)], we use the galaxy number 
counts, g(< m) = ioO-33('«-24)-2.32 arcsec-2 (Hogg et al 1997; 
iBeckwith et al.l 12006). For the XRT error circle we infer 
Pec ~ 0.1, while for the more precise JVLA and Chandra posi- 
tions we find P^c ~ 0.006 and w 0.001, respectively. Thus, the 
sub-arcsecond radio and X-ray positions enable a robust asso- 
ciation with the galaxy, while the XRT position alone would 
have led to marginal confidence. This is essential since the 
detection of the host galaxy in the F606W filter limits the red- 
shift to z < 4, indicating that GRB 110709B is dark due to 
extinction and not a z > 18 origin. 

At z ^ 2, the observed F606W band samples the rest-frame 
UV emission from the host galaxy. Based on the observed 
flux density we in fer a star formation rate of SFR 1 Mq yr"' 
(iKennicuttI 1 998h . not corrected for possible galactic-scale ex- 
tinction. 

2.6.2. GRB 11 121 5A 

We obtained z-band observations of GRB 111215A with 
GMOS on the Gemini-North 8-m telescope on 2011 De- 
cember 23.21 UT (St w 7.6 d), with a total exposure time 



of 2160 s in 0.9" seeing. Astrometric matching to the 
2MASS catalog using 30 common sources reveals a galaxy 
within the enhanced XRT eiTor circle at RA=23''18™13.317', 
Dec=+32°29'38.77", with an uncertainty of about 0.16" in 
each coordinate (absolute) and an uncertainty of about 0.06" 
in the galaxy centroid; see Figure|5] We also obtained imaging 
observ ations with the Low-Resolution Imaging Spectrometer 
(LRIS; IOke et al.ll995 ) on the Keck 1 10-m telescope on 2012 
July 15.20 UT (St = 212.6 d) in the g- and /-band filters. 

Taking into account the uncertainty in the JVLA radio after- 
glow position (0.06" in each coordinate), the resulting over- 
all uncertainty in the location of the afterglow on the Gemini 
image is 0.26" radius (la). Similarly, the uncertainty in the 
CARMA mm afterglow position (0.08" in each coordinate) 
leads to an overall uncertainty of 0.27" radius (la). The re- 
sulting offset between the host centroid and radio positions is 
0.30 ±0.26" 

Photometry of the galaxy relative to the SDSS catalog re- 
sults in a brightness of nig = 24.50 ±0.15 mag and m,- = 
23 .64 ±0.05 mag corrected for Galactic extinction. The prob- 
ability of chance coincidence for this galaxy within the XRT 
error circle is fee ~ 0.01, with a slightly lower probability of 
w 7 X 10""^ using the radio position, dominated by the phys- 
ical extent of the galaxy of about 1". Given the low chance 
coincidence probability we consider this galaxy to be the host 
ofGRB 111215A. 

To determine the redshift of the host we obtained an 1800 s 
specti-um with LRIS, using a 0.7" slit with the 400/3400 grism 
in the blue arm and the 400/8500 grating in the red arm, pro- 
viding effectively continuous wavelength coverage from the 
atmospheric cutoff to 10280 A (with a total throughput within 
50% of peak to about 10030 A). The slit was oriented to 
cover the host galaxy, as well as the fainter extended source to 
the northeast. The spectra were reduced using standard tech- 
niques implemented in a custom pipeline. 

We detect continuum emission from both galaxies on the 
blue side, and from the host galaxy on the red side. We 
identify only one emission feature from the nearby galaxy, a 
marginally resolved line centered at 8019 A. Interpreted as the 
[O II]A3727 doublet, this indicates a redshift'^ of z = 1.152. 
We do not identify any emission lines from the host galaxy. 
Given that the galaxy is reasonably bright in g- and /-band 
and is expected to be actively star-forming, the lack of emis- 
sion lines indicates that z > 1.8 (from the non-detection of 
[O II]A3727 to 10030 A) and z < 2.7 (from the absence of a 
break due to the Lyman limit to 3400 A). 

At z ^ 1.8-2.7, the observed /- and g-band fluxes trace 
the host galaxy rest-frame UV emission. The observed color 
of g-i ^ 0.85 ± 0.20 mag is indicative of extinction; for ex- 
ample, it is well-matched to the observed color of Arp 220, a 
local ULIRG. Using the observed /-band flux density, which 
is less susceptible to extinction corrections, we infer SFR w 
15-30 M© yr"' (z = 1.8-27), which should be considered as 
a minimum value due to the uncertain extinction. 

3. AFTERGLOW MODELING AND REST-FRAME EXTINCTION 

GRBs 110709B and 111215A are dark in both the optical 
and near-IR bands due to rest-frame dust extinction; a high- 
redshift origin is ruled out by the association with optically- 
detected host galaxies. It is also of note that for extreme red- 

We also examined the possibility of alternative line identifications and 
redshifts, but these would imply the detection of other spectral lines else- 
where that are not observed. 
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shifts (z > 18), the A^H.int impHed by the observed excess in 
the X-ray spectra would be very high (e.g. > 10^'* cm"^), 
further supporting the rejection of an extreme redshift for ei- 
ther event. In addition to providing accurate localizations and 
hence secure host associations, the radio data also allow us to 
determine the properties of the bursts and their local environ- 
ments. This not only provides a more robust measure of the 
required rest-frame extinction than using /3ox alone (which 
requires an assumption about the location of the synchrotron 
cooling frequency), but it also allows us to compare the ex- 
plosion properties of dark and optically-bright bursts. 

We model the radio and X-ray data for GRBs 110709B 
and 1 1 1215 A using the standard afterglow synchrotron model 
jGranot & Sai'i 2002; Sari et al. 1999) to: (i) determine the 
expected optical and near-lR brightness, and hence the re- 
quired level of rest-frame extinction given the observed lim- 
its; and (ii) determine the properties of the bursts and their 
circumburst environments. We follow the standard assump- 
tions of synchrotron emission from a power-law distribution 
of electrons (A^(7) oc 7"^ for 7 > 7,,,) with constant fractions 
of the post-shock energy density imparted to the electrons (eg) 
and magnetic fields (e^). The additional free parameters of 
the model are the isotropic-equivalent blast-wave kinetic en- 
ergy (Eksso), the circumburst density, (parametrized as n for 
a constant density medium: ISM; or as A for a wind medium 
with p(r) =Ar~^), and a jet break time (f^). To determine the 
opening angle (O j), we use the conversions from tj given by 
ISariet al.1 (Il999l) and lChevaUer & Lil (|2000|) . with the appro- 
priate dependence on iiK.iso and the circumburst density (n or 
A, respectively). 

Using the time evolution of the synchrotron spectrum for 
both the ISM and wind density profi l es in the pre- and pos t-jet 
break phases (IChevaUer & LillioOOl; iGranot & Sarill2cio2l) we 
simultaneously fit all X-ray and radio observations for each 
burst. The resulting best-fit parameters are summarized in Ta- 
bles|4]and|5]usingredshifts of z = 1,2,3,4 for ORB 110709B 
and z = 2 for GRB 111215A. In both cases we find that the 
most stringent constraint on the rest-frame extinction are pro- 
vided by the observed /T-band limits. We use the SMC extinc- 
tion curve to determine the rest-frame extinction, although the 
Milky Way extinction curve gives similar results. 

For GRB 110709B we find that the data favor a wind en- 
vironment, with a jet break at tj w 3.1-3.6 d, correspond- 
ing to an opening angle, 9j, w 9-23°; the jet break oc- 
curs earlier and the resulting opening angle is narrower as 
the redshift increases from z = 1 to z = 4. The resulting 
beaming-corrected energies are Ek ~ (0.4- 1.7) x 10^' erg 
and w (5.4- 10.6) x 10^' erg (in the observed 20-5000 
keV fluence of the initial event). The circumburst density is 
characterized by a mass loss rate of M w (4- 8) x 10"^ Mq 
yr"' for a wind velocity of v„, = lO-' km s"' (i.e.. A* ss 4- 8). 
Finally, using the near-IR /Q-band limit we find that the re- 
quired extinction ranges from A{j,°^' > 10.5 mag at z = 1 to 
> 3.5 mag at z = 4. We note that the best-fit ISM model also 
requires a high density of 10- 100 cm"-' at z 1-4. 

For GRB 1 1 12 15 A we again find that a wind environment 
provides a better fit, with a jet break at tj « 12 d, cor- 
responding to an opening angle of 6j w 24°. The result- 
ing beaming-corrected energies are Zs^: ~ 3.7 x 10^' erg and 
E^ w 3.9 X 10^' erg (in the observed 15-150 keV band). The 
circumburst density is characterized by a substantial mass loss 
rate of M « 2 x 10""* M0 yr"' (for a wind velocity of v„ =10^ 
km s"'). Finally, using the near-lR /Tj-band limit we find a re- 



quired extinction of Ay°^' ^8.5 mag. We note that the best-fit 
ISM model also requires a density of ^ 100-300 cm"^. 
4. DISCUSSION 
4.1. Extinction and Neutral Hydrogen Column Density 

The large rest-frame extinction that is required to explain 
the lack of optical and near-IR emission from GRBs 1 10709B 
and 111215A is uncommon amongst known GRBs. In Fig- 
ure |6] we plot the minimum values of A{^°^' for the two 
bursts as a function of redshift, along with several compar- 
ison samples from the literature of optically-bright bursts 
and previous dark bursts. In general, optically-bright GRBs 
have inferred extinction values of < 1 ma g, with a mean 
of ^ 0.2- .3 mag (Kannetal. 2006; Sc hadv et al.l 120071; 
lPerley etal.ll200 9b; Kruhleretal. 2011) . Dark GRBs spa n 
Af'' « 1-6 mag (Pe rley et al 2009b; .Kriihler et al.ll20lTl) . 
with the largest extinction values of ~ 5-6 mag inferred 
for GRBs 061222, 070306, and 070521 (iPerley et al.ll2009bl: 
Kruhleretal. 2011). In this context, the extinction values 
measured here are at the top of the distribution, Ay°'*' ^5.3 
mag ( 1 1 0709B) and > 8 .5 mag ( 1 1 1 2 1 5 A) at z = 2. 

In Figure |7] we plot the intrinsic neutral hydrogen col- 
umn densities, A^n.int, from our X-ray analysis as a func- 
tion of redshift, in comparison with the population of all 
Swift long GR Bs w i th kn own redshifts (up to December 
2010; Margu tti etaH I20TI) . We find that GRBs 110709B 
and 1 1 1215A lie at the upper end of the distribution for long 
GRBs. In particular, at a fiducial redshift of z = 2 they have 
log(A^H.int) ~ 22.1 and w 22.5, respectively, in comparison 
to the median value of all detections and upper limits of 
(log(A^H.mt)) ~ 21.7. This result confirms recent suggestions 
that dark GRBs have systematically larg er neutral hydrogen 
columns, generally w ith log(A^H,int) ^ 22 (iPerley et al.ll2009al; 
ICampana et al.ll2012l) ; see Figure |7] Indeed, only a handful of 
events in the Swift sample have values of A^e.int that are com- 
parable to that of GRB 1 1 12 15 A. The same conclusion is true 
for GRB 1 10709B if it resides at the upper end of the allowed 
redshift distribution (z 3 - 4), although at z ~ 1 its neutral 
hydrogen column is typical of the overall long GRB popula- 
tion. 

A correlation between extinction and Nh is known to ex- 
ist in the Milky Wa y and the Magellanic Clo u ds, with Nh ~ 
2 X 1 0^' Ay cm"2 (jPredehl & Schiiiita [T99l iGuver & OzeJ 
l2009l) . As shown in Figure H] optically-bright GRBs gener- 
ally have lower values of A5)°^' than would be inferred from 
this relation (e.g., Galama & Wiiers 2001), with mean values 
of A^aint « 5 X 10^1 cm"2 and A^"'' « 0.2-0.3 mag; the ex- 
pected extinction based on the Galactic relation is A^°^' « 2.5 
mag. Dark GRBs also generally have lower extinction than 
expected, with A^Hant « (5-50) x 10^' cm"^ and A^°'*' « 1 - 5 
mag, whereas extinctions of A^°^' « 5 - 25 mag would be ex- 
pected. 

For the bursts presented here we find potentially different 
results. GRB 111215A agrees with the Galactic relation if it 
resides at the lower redshift bound (z ~ 1.8), in which case 
A^H.int ~ 2.5 X 10^2 cm"2 and A^°'' > 11 mag. At the higher 
redshift bound (z ~ 2.7) we find A^H.int ~ 5 x 10^^ cm"^ and 
^host > jj^^g (yyji-jj expected Ay°'*' « 23 mag). This is 
still consistent with the Galactic relation since we only place 
a lower bound on Aj)?". Similarly, GRB 1 10709B crosses the 
Galactic relation at z ~ 2, with A^H.int ~ 1.2 x 10^^ cm"^ and 
Ay°^' > 5.3 mag. However, at z < 2 it has larger extinction than 
expected; for example, at z = 1 we measure Ay°'*' > 10.5 mag 
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compared to the expected value of « 2.5 mag. At z > 2 the 
situation is similar to GRB 11 12 15 A, namely the minimum 
value of Ay°'*' is below the expected Galactic relation, but this 
is only a lower limit. We therefore conclude that the dark 
bursts with the largest known extinction values are potentially 
in line with the Galactic Nn-Ay relation, although they may 
also reside below the expected correlation in line with the bulk 
of the GRB population. 

4.2. Explosion and Circumburst Properties of Dark Bursts 

A few previous dark GRBs have been localized to sub- 
arcsecond precision in the radio. GRB 970828 was detected 
in a single radio observation with /v(8.46GHz) = 147 ± 33 
/iJy (4.5(t), but was r iot detected in observa tions only 1 d 
earlier and 2 d later jDiorgovski et alj |2001|) . This single 
marginal detection was used to refine the position from the 
X-ray region of 10" radius, and to claim an association with a 
galaxy at z = 0.958. The inferred extinction at this redshift is 
A^?" > 3.8 mag. GRB 000210 was localized to sub-arcsecond 
precision with Chandra, and was also marginally detected in 
the single epoch in the rad io with /v (8. 46 GHz) = 93 ± 21 
/xJy (4.4(t; iPiro et aLll2002h . The associated host galaxy has 
z = 0.846, and the inferred extinction is Ajj,™' « 0.9-3.2 mag. 
GRB 0208 19 was detected at high significance in several 
epochs of radio observations, leading t o an association wit h 
a face-on spiral galaxy at z = 0.410 (Jakobsson et alj|2005h . 
The required extinction was a mild Ay°^' w 0.6- 1 .5 mag. Fi- 
nally, GRB 051022 was localized using radio and milUmeter 
observations to a host gala xy at z = 0.809, requiring an extinc- 
tion o f Ajf > 3.5 mag ( Castro-Tirado et al. 20ol iRol et alj 
l2007h . Broad-band modeling indicates a beaming-corrected 
energy of ~ 10^' erg, and a circumburst density characterized 
by M w 3 X 10"^ Mq yr'^ (iRol et al.ll2007h . 

The beaming-corrected energy inferred for GRBs 1 10709B 
and 111215A, E^+Ek w (7-9) x 10^' erg, is not unusual 
for long GRBs, although it is some what higher than the 
median value of w 3 x 1 0^^ erg (e.g.. iPanaitescu & Kumaii 
I2002t iBerger et al.ll2003bl: ICenko et al.ll201 lb . On the other 
hand, the inferred densities are substantial, with inferred 
mass loss rates of about 6 x 10"^ M© yr"' (110709B) and 
2 x 10""* M© yr"' (111215A). At the radii appropriate for 
the jet break times of the two bursts, r « (1.2- 1.4) x 10'^ 
cm (C hevalier & Lill2000l) . these mass loss rates correspond 
to particle densities of about 100 and 350 cm"^, respectively, 
higher than typical va lues for long GRBs , which are ^ 0. 1-3 
cm"^ (e.g.. IChevalier & Li 2000; Panaitescu & Kumadl2002h . 
If confirmed with future detailed observations of dark bursts, 
the high densities may be indicative of larger mass loss rates 
for long GRBs in dusty, and hence metal rich environments. 
Such a trend may be expected if the mass loss is driven by 
radiation pressure mediated by metal lines dVink & de Koterl 
l2005h . 

4.3. Burst Durations 

The durations of GRB110709B (r^ w 850 s) and 
GRB 111215A (Tgo « 800 s) are extremely long. To assess 
whether this type of unusually long duration correlates with 
optical/near-IR darkness, we collect all Swift bursts with du- 
rations of > 500 s, excluding nearby sub-energetic GRBs 
(e.g., XRF060218, GRB 100316D). The sample includes 9 
events'"* of which three bursts ai-e dark (060929, 110709B, 

'* These are GRBs 041219A, 050820A, 060123, 060124, 060929, 091024, 
110709B, 111016A, and 111215A. 



11 12 15 A), five have optical and/or near-IR detections with 
no clear evidence for extinction, and one event lacks rapid 
follow-up observations. Thus, at least ~ 1/3 of these events 
are dark, similar to the fraction of dark bursts in the overall 
long GRB population. We therefore conclude that there is 
no obvious correlation between unusually long duration and 
darkness, with the caveat that the current sample size is small. 

5. CONCLUSIONS 

Using X-ray, optical/near-IR, and radio observations we 
have demonstrated that: (i) GRBs 1 10709B and 1 1 12 15 A are 
dark bursts; (ii) they are robustly associated with galaxies at 
z < 4 (110709B) and z « 1.8-2.7 (111215A); (iii) they re- 
quire unusually large rest-frame extinction of Ajf**' > 5.3 mag 
(1 10709B) and > 8.5 mag (1 1 1215A) at z = 2; (iv) they ex- 
hibit commensurately large neutral hydrogen column densi- 
ties in their X-ray spectra, A^H.int ~ (1 - 3) x 10^^ cm"^, which 
at z ^ 2 are consistent with the Galactic Nh - Ay relation, un- 
like the overall long GRB population; and (v) their circum- 
burst environments on a sub-parsec scale are shaped by large 
progenitor mass loss rates of w (6-20) x 10"^ M© yr"'. 

Radio observations played a critical role in this study for 
three reasons. First, they provided sub-arcsecond positions, 
which led to a secure identification of the host galaxies. This 
allowed us to distinguish the extinction scenario from a high- 
redshift origin. Second, the combination of radio and X-ray 
data allowed us to robustly determine the required extinction, 
instead of simply assuming an optical to X-ray spectral index. 
Indeed, in cases with only X-ray data the unknown location 
of the synchrotron cooling frequency prevents a unique de- 
termination of the extinction. Finally, the radio and X-ray 
data allowed us to determine the burst parameters, including 
the geometry, beaming-corrected energy, and the circumburst 
density. We find that the energy scale for GRBs 110709B 
and 1 1 12 15 A is similar to the overall population of optically- 
bright GRBs. However, the inferred mass loss rates are larger 
by about an order magnitude compared to optically-bright 
bursts, potentially indicating that GRB progenitors in dusty 
environments have stronger metal line driven winds. 

The increased sensitivity of the JVLA played an important 
role in the study and utilization of the radio afterglow. ALMA 
will provide similar capabilities, with an expected spatial res- 
olution of about 0.04-5" at 100 GHz (depending on con- 
figuration), and hence an assured sub-arcsecond centroiding 
accuracy for signal-to-noise ratios of > 10 even in the most 
compact configuration. As demonstrated here, Chandra ob- 
servations can also provide sub-arcsecond positions, but such 
observations may not be available for all dark GRB candidates 
(e.g., GRB 111215A). In addition, Chandra data will not en- 
hance the ability to determine the required extinction through 
broad-band modeling. 

The JVLA and ALMA will also allow us to study the host 
galaxies of GRBs 1 10709B and 1 1 12 15 A (as well as those of 
past and future dark bursts) in greater detail than optical/near- 
IR studies alone. In particular, they will address the ques- 
tion of highly-obscured star formation in these galaxies, and 
may even lead to redshift determinations through the detec- 
tion of molecular lines. Existing observations with the VLA 
and the MAMBO and SCUBA bolometers (in the small sam- 
ples published to-date) have led to only a few detections of 
GRB hosts (Berger et al. 200 3^, with no obvious preference 
for dark burst hosts (Bar nard et al.ll20()3h . However, the en- 
hanced sensitivity of ALMA will allow for detections even 
with star formation rates of tens of M© yr"' at z ^ 2. The com- 
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bination of detailed host galaxy properties from rest-frame ra- 
dio to UV, coupled with measurements of the local environ- 
ments of dark bursts through broad-band afterglow modeUng 
will shed light on the location of the obscuring dust (interstel- 
lar vs. circumstellar) and the impact of metallicity on GRB 
progenitor formation. 
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TABLE 1 

GROND Observations of GRB 1 10709B 



UTDate 


St 


Timespan 




Upper Limits 




Timespan 


Upper Limits 




(d) 


(s) 








(s) 










g' 


r' i' 


z' 


J H Ks 


201 1 Jul 9.949 


0.051 


585 


20.76 


20.53 20.32 


20.49 


621 


20.03 19.02 18.79 


201 1 Jul 10.009 


0.111 


2311 


24.26 


24.57 23.83 


23.61 


2335 


20.83 20.00 19.27 


201 1 Jul 10.045 


0.147 


1633 


22.61 


22.80 22.39 


21.12 


1640 


19.83 18.83 18.66 


2011 Jul 10.003" 


0.105 


5275 


23.63 


23.89 23.28 


23.09 


5282 


21.16 19.83 19.40 



Note. — Date specifies the mid-time of the observations, and timespan is the duration of the observations. AH 
magnitudes are upper limits in the AB system, uncorrected for Galactic extinction. 
" Results from stacking all three epochs. 
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TABLE 2 

JVLA Observations of GRB 1 10709B 



UT Date 


St 


V 






(d) 


(GHz) 




2011Jul 11.97 


2.07 


5.8 


190 ±14 


2011 Jul 12.92 


3.02 


5.8 


250 ± 27 


2011 Jul 16.88 


6.98 


5.8 


310±18 


2011 Jul 21.97 


12.07 


5.8 


210d=14 


2011 Jul 29.96 


20.06 


5.8 


98± 17 


2011 Aug 19.85 


40.59 


5.8 


98±6 


2011 Sep 17.70 


69.72 


5.8 


< 38" 


2011 Jul 22.01 


12.11 


21.8 


< 135 



Note. — Dates specify the mid-time of the 
observations. 
" Upper limits are 3a. 
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TABLE 3 

Radio Observations of GRB 111215A 



Telescope 


UT Date 


5t 


V 








(d) 


(GHz) 




JVLA 


201 1 Dec 17.00 


1.41 


4.9 


< 115" 








6.7 


< 93 




zuii uec zj.uj 


/ .40 


A Q 


^ / U It zo 








6.7 


460 ± 34 




201 1 Dec 29.03 


13.44 


4.9 


300 ±21 








6.7 


540 lb 19 




701 1 ^ 1 on 


1 S A1 




890 -1- '^O 








6. / 


1 1 OA _L ■lA 

1 1 oO ± 30 




2012 Jan lu.yy 


20.40 


4.9 


CAA J_ I'l 

jyO ± 33 








0./ 


OOA _L 

o20 ± 3d 




iKjvz Jan 2y.o4 


43.2!) 


A O 

4.9 


/f^A _L CO 

420 ± j2 








0. / 


/lAA J_ CC 

4U0 iL jD 




2012 Mar 12. o2 


00.23 


4.9 


330 ± 39 








0./ 


■2AA 1 "> A 

300 ± 34 




2011 Dec 31.00 


15.41 


8.4 


1340 ±63 




2012 Jan 10.99 


26.40 


8.4 


880 ± 120 




2012 Jan 29. o4 


45.25 


8.4 


470di 81 




2012 Mar 12. o2 


88.23 


8.4 


1 /I A _1_ AC 

240 ± 95 




zuii uec io.y^ 


^ ^< 

J.JO 


1 Q 1 










24.4 


1 990 ± 47 




2011 Dec 27.0iS 


11.49 


19.1 


lo/O ± 47 








24.4 


1650 di 74 




OA1 1 r\,^ ,-. O 1 AA 

2011 Dec 31.00 


15.41 


19.1 


1 CIA _1_ A A 

1530 ± 49 








24.4 


1370 ±62 




9019 Tnn 10 QQ 


9fi ,10 


10 1 










24.4 


730 ±78 




2012 Jan 29.84 


45.25 


19.1 


500 ±66 








24.4 


460 ±69 




2012 Mar 12.82 


88.23 


19.1 


214 ±66 








24.4 


180 ±72 


CARMA 


2011 Dec 17.05 


1.46 


93 


3400 ±290 




2011 Dec 20.02 


4.43 


93 


2500 ±330 




2011 Dec 22.04 


6.45 


93 


2500 ±300 




2011 Dec 24.10 


8.51 


93 


2200 ±270 




2011 Dec 26.01 


10.42 


93 


2200 ± 270 




2012 Jan 1.14 


16.55 


93 


1500 ±340 


SMA 


2011 Dec 18.30 


2.71 


230 


<2600 



Note. — Dates specify the mid-time of the observations. 
" Upper limits are 3it. 
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TABLE 4 

Results of Broad-band Afterglow Modeling 

1 OR GRB 1 I 07()9B 



Parameter 


z=l 


1 = 2 


2 = 3 


z = 4 


£'K,iso,52 (erg) 


0.5 


2.8 


7.1 


13.2 


A* (5 X 10" g cm- 


1) 4.0 


5.7 


7.2 


8.5 




0.18 


0.11 


0.08 


0.07 




0.002 


0.002 


0.002 


0.002 


P 


2.07 


2.06 


2.05 


2.05 


tj (d) 


3.65 


3.45 


3.27 


3.11 


9j{rad) 


0.40 


0.25 


0.19 


0.16 


Ek,51 (erg) 


0.4 


0.9 


1.3 


1.7 


E-fAsoSl" (erg) 


7.0 


26 


52 


83 


£■^,^51 (erg) 


5.4 


8.0 


9.4 


10.6 


Ai/ * (mag) 


> 10.5 


>5.3 


>4.4 


?.3.4 



Note. — See §3 for description of model parameters. 
" Using tlie Konus-WIND 20-5000 keV fluence. 
* Using the SMC extinction ciu-ve. 
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TABLE 5 

Results of Broad-band Afterglow 

MODliLiNG l-OR GRB I 1 I 2 I 5 A 



Parameter 


z = 2 


z = 3 


£^K,iso,52 (erg) 


4.3 


10.6 


A* (5 X 10" g cm-i) 


19 


26 


£« 


0.22 


0.18 


es 


2.2 X 10-4 


1.5 X 10-4 


P 


2.35 


2.35 


t] (d) 


12 


12 


6», (rad) 


0.42 


0.33 


£k.51 (erg) 


3.7 


5.7 


£■7,150,52" (erg) 


4.5 


9.0 


£7,51 (erg) 


3.9 


4.9 


Av* (mag) 


>8.5 


>6.8 



Note. — See §3 for description of model pa- 
rameters. 

" Using the fluence in the SwiftlBNY 15 - 150 keV 
band. 

* Using the SMC extinction curve. 
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Time since burst (d) 



Fig. 1 . — X-ray (black) and radio (blue) light curves of GRB 1 10709B, with near-IR limits (red triangles). For the purpose of display, the XRT data are binned 
by orbit at 5f < 2 d and by multiple orbits thereafter. The Chandra 3cr upper limit at (5f Ri 1 13 d is indicated with a black arrow. Afterglow model fits in each band 
are shown with soHd lines. To satisfy the near-IR upper limits we require Ay > 3.4-10.5 mag (for z = 4 to 1; dashed red line). The derived afterglow parameters 
for GRB 1 10709B are listed in Table|4] 
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Fig. 2. — X-ray (black) and radio (blue) light curves of GRB 1 1 1215A, with near-IR limits (red triangle). Afterglow model fits in each band are shown with 
solid lines. To explain the near-IR upper limit we require Ay > 8.5 mag (dashed line, assuming z = 2). The derived afterglow parameters are listed in Table[5] 
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Gemini-North / GMOS / i-band 
2011 December 23.21 UT 



I 



EVLA/ 
CARMA 





XRT ,r 



2" 




Fig. 5. — Gemini/GMOS /-band image of the host galaxy of GRB 1 1 1215A. The XRT error circle (1.4" radius; 90% containment) is indicated with the black 
circle. The JVLA and CARMA afterglow positions are shown with a blue circle (0.13"; lo" systematic) coinciding with a galaxy that we consider to be the host 
of GRB 111215A. 
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' Dark Bursts - solid symbols ' 
Other GRBs - open symbols 
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Redshift 
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GRB 110709B 
GRB111215A 
GRB050401 
Kruhler et al. 2011 
Perley/Cenko 



□ GRB 090423 

O Wijers & Galama 2001 

O Schady et al. 2007 

O Kann et al. 2006 

o Kann et al. 2010 

O Perley/Cenko 




Fig. 6. — Rest-frame extinction (A^ versus redshift for GRBs 110709B (blue triangles) and 111215A (red triangle); for both bursts these are lower limits. 
Previous dark bursts are marked with solid symbols. Triangles indicate limits on AJ,"*' and/or redshift. Allowed values for GRBs without firm redshifts are 
indicated by black lines jPerlev et aU20093) . Optically-bright GRBs are marked with open symbols; these bursts gener ally have AJ,"'" < 1 mag. GRBs 110709B 
and 1 1 1215A exhibit some of the highest extinction values to date. 




2 4 6 8 

Redshift 

Fig. 7. — Intrinsic neutral liydrogen column density (A'n.int) inferred from X-ray observations versus redshift for GRBs 110709B and 111215A (contours). 
The allowed redshift ranges for each burst are marked by the gray shaded regions. For comparison, we overlay data for long GRBs in the Swift sample up to 
December 2010 ( Mar gutti et al . 2012); dark bursts are marked by blue/black circles, while optically-bright events are marked by gray circles. GRBs 110709B 
and 1 1 1215A clearly lie in the upper portion of the A^Hjint distribution of long GRBs, as expected for dark bursts. 
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Kann et al. 2006 
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Kann et al. 2010 
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Perley et al. 2009 



10 20 -2 10" 

N„ . ^ (10^^ cm ^) 
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Fig. 8. — Rest-frame extinction (AJ,"*') versus intrinsic neutral liydrogen column density (Wnant) inferred from X-ray observations for GRBs 110709B (blue 
triangles) and 1 1 1215A (red triangle). For compaiison, we plot the same sample of bursts shown in Figure[6] The diagonal Unes indicate the empirical Ay —Nh 
relations for the Milky Way and Magellanic Clouds. Most GRBs exhibit lower values of Ay""^ than expected from the local relations, but GRBs 1 10709B and 
1 1 1215A may be consistent with these relations. 



